Specific oncogene activation of the cell of origin in mucosal melanoma.

Mucosal melanoma (MM) is a deadly cancer derived from mucosal melanocytes. To test the consequences of MM genetics, we developed a zebrafish model in which all melanocytes experienced CCND1 expression and loss of PTEN and TP53. Surprisingly, melanoma only developed from melanocytes lining internal organs, analogous to the location of patient MM. We found that zebrafish MMs had a unique chromatin landscape from cutaneous melanoma. Internal melanocytes could be labeled using a MM-specific transcriptional enhancer. Normal zebrafish internal melanocytes shared a gene expression signature with MMs. Patient and zebrafish MMs have increased migratory neural crest gene and decreased antigen presentation gene expression, consistent with the increased metastatic behavior and decreased immunotherapy sensitivity of MM. Our work suggests the cell state of the originating melanocyte influences the behavior of derived melanomas. Our animal model phenotypically and transcriptionally mimics patient tumors, allowing this model to be used for MM therapeutic discovery.

Single-cell RNA sequencing data imputation using bi-level feature propagation.

Single-cell RNA sequencing (scRNA-seq) enables the exploration of cellular heterogeneity by analyzing gene expression profiles in complex tissues. However, scRNA-seq data often suffer from technical noise, dropout events and sparsity, hindering downstream analyses. Although existing works attempt to mitigate these issues by utilizing graph structures for data denoising, they involve the risk of propagating noise and fall short of fully leveraging the inherent data relationships, relying mainly on one of cell-cell or gene-gene associations and graphs constructed by initial noisy data. To this end, this study presents single-cell bilevel feature propagation (scBFP), two-step graph-based feature propagation method. It initially imputes zero values using non-zero values, ensuring that the imputation process does not affect the non-zero values due to dropout. Subsequently, it denoises the entire dataset by leveraging gene-gene and cell-cell relationships in the respective steps. Extensive experimental results on scRNA-seq data demonstrate the effectiveness of scBFP in various downstream tasks, uncovering valuable biological insights.

Biophysical characterization of high-confidence, small human proteins.

Significant efforts have been made to characterize the biophysical properties of proteins. Small proteins have received less attention because their annotation has historically been less reliable. However, recent improvements in sequencing, proteomics, and bioinformatics techniques have led to the high-confidence annotation of small open reading frames (smORFs) that encode for functional proteins, producing smORF-encoded proteins (SEPs). SEPs have been found to perform critical functions in several species, including humans. While significant efforts have been made to annotate SEPs, less attention has been given to the biophysical properties of these proteins. We characterized the distributions of predicted and curated biophysical properties, including sequence composition, structure, localization, function, and disease association of a conservative list of previously identified human SEPs. We found significant differences between SEPs and both larger proteins and control sets. Additionally, we provide an example of how our characterization of biophysical properties can contribute to distinguishing protein-coding smORFs from non-coding ones in otherwise ambiguous cases.

Single-cell dissection of the human motor and prefrontal cortices in ALS and FTLD.

Amyotrophic lateral sclerosis (ALS) and frontotemporal lobar degeneration (FTLD) share many clinical, pathological, and genetic features, but a detailed understanding of their associated transcriptional alterations across vulnerable cortical cell types is lacking. Here, we report a high-resolution, comparative single-cell molecular atlas of the human primary motor and dorsolateral prefrontal cortices and their transcriptional alterations in sporadic and familial ALS and FTLD. By integrating transcriptional and genetic information, we identify known and previously unidentified vulnerable populations in cortical layer 5 and show that ALS- and FTLD-implicated motor and spindle neurons possess a virtually indistinguishable molecular identity. We implicate potential disease mechanisms affecting these cell types as well as non-neuronal drivers of pathogenesis. Finally, we show that neuron loss in cortical layer 5 tracks more closely with transcriptional identity rather than cellular morphology and extends beyond previously reported vulnerable cell types.

FinaleMe: Predicting DNA methylation by the fragmentation patterns of plasma cell-free DNA.

Analysis of DNA methylation in cell-free DNA reveals clinically relevant biomarkers but requires specialized protocols such as whole-genome bisulfite sequencing. Meanwhile, millions of cell-free DNA samples are being profiled by whole-genome sequencing. Here, we develop FinaleMe, a non-homogeneous Hidden Markov Model, to predict DNA methylation of cell-free DNA and, therefore, tissues-of-origin, directly from plasma whole-genome sequencing. We validate the performance with 80 pairs of deep and shallow-coverage whole-genome sequencing and whole-genome bisulfite sequencing data.

Gamma entrainment using audiovisual stimuli alleviates chemobrain pathology and cognitive impairment induced by chemotherapy in mice.

Patients with cancer undergoing chemotherapy frequently experience a neurological condition known as chemotherapy-related cognitive impairment, or "chemobrain," which can persist for the remainder of their lives. Despite the growing prevalence of chemobrain, both its underlying mechanisms and treatment strategies remain poorly understood. Recent findings suggest that chemobrain shares several characteristics with neurodegenerative diseases, including chronic neuroinflammation, DNA damage, and synaptic loss. We investigated whether a noninvasive sensory stimulation treatment we term gamma entrainment using sensory stimuli (GENUS), which has been shown to alleviate aberrant immune and synaptic pathologies in mouse models of neurodegeneration, could also mitigate chemobrain phenotypes in mice administered a chemotherapeutic drug. When administered concurrently with the chemotherapeutic agent cisplatin, GENUS alleviated cisplatin-induced brain pathology, promoted oligodendrocyte survival, and improved cognitive function in a mouse model of chemobrain. These effects persisted for up to 105 days after GENUS treatment, suggesting the potential for long-lasting benefits. However, when administered to mice 90 days after chemotherapy, GENUS treatment only provided limited benefits, indicating that it was most effective when used to prevent the progression of chemobrain pathology. Furthermore, we demonstrated that the effects of GENUS in mice were not limited to cisplatin-induced chemobrain but also extended to methotrexate-induced chemobrain. Collectively, these findings suggest that GENUS may represent a versatile approach for treating chemobrain induced by different chemotherapy agents.

New approaches for understanding the potential role of microbes in Alzheimer's disease.

Alzheimer's disease (AD) involves a complex pathological process that evolves over years, and its etiology is understood as a classic example of gene-environment interaction. The notion that exposure to microbial organisms may play some role in AD pathology has been proposed and debated for decades. New evidence from model organisms and -omic studies, as well as epidemiological data from the recent COVID-19 pandemic and widespread use of vaccines, offers new insights into the "germ hypothesis" of AD. To review new evidence and identify key research questions, the Duke/University of North Carolina (Duke/UNC) Alzheimer's Disease Research Center hosted a virtual symposium and workshop: "New Approaches for Understanding the Potential Role of Microbes in Alzheimer's disease." Discussion centered around the antimicrobial protection hypothesis of amyloid accumulation, and other mechanisms by which microbes could influence AD pathology including immune cell activation, changes in blood-brain barrier, or direct neurotoxicity. This summary of proceedings reviews the content presented in the symposium and provides a summary of major topics and key questions discussed in the workshop.

Disrupting CD38-driven T cell dysfunction restores sensitivity to cancer immunotherapy.

A central problem in cancer immunotherapy with immune checkpoint blockade (ICB) is the development of resistance, which affects 50% of patients with metastatic melanoma1,2. T cell exhaustion, resulting from chronic antigen exposure in the tumour microenvironment, is a major driver of ICB resistance3. Here, we show that CD38, an ecto-enzyme involved in nicotinamide adenine dinucleotide (NAD+) catabolism, is highly expressed in exhausted CD8+ T cells in melanoma and is associated with ICB resistance. Tumour-derived CD38hiCD8+ T cells are dysfunctional, characterised by impaired proliferative capacity, effector function, and dysregulated mitochondrial bioenergetics. Genetic and pharmacological blockade of CD38 in murine and patient-derived organotypic tumour models (MDOTS/PDOTS) enhanced tumour immunity and overcame ICB resistance. Mechanistically, disrupting CD38 activity in T cells restored cellular NAD+ pools, improved mitochondrial function, increased proliferation, augmented effector function, and restored ICB sensitivity. Taken together, these data demonstrate a role for the CD38-NAD+ axis in promoting T cell exhaustion and ICB resistance, and establish the efficacy of CD38 directed therapeutic strategies to overcome ICB resistance using clinically relevant, patient-derived 3D tumour models.

Rabies virus-based barcoded neuroanatomy resolved by single-cell RNA and in situ sequencing.

Mapping the connectivity of diverse neuronal types provides the foundation for understanding the structure and function of neural circuits. High-throughput and low-cost neuroanatomical techniques based on RNA barcode sequencing have the potential to map circuits at cellular resolution and a brain-wide scale, but existing Sindbis virus-based techniques can only map long-range projections using anterograde tracing approaches. Rabies virus can complement anterograde tracing approaches by enabling either retrograde labeling of projection neurons or monosynaptic tracing of direct inputs to genetically targeted postsynaptic neurons. However, barcoded rabies virus has so far been only used to map non-neuronal cellular interactions in vivo and synaptic connectivity of cultured neurons. Here we combine barcoded rabies virus with single-cell and in situ sequencing to perform retrograde labeling and transsynaptic labeling in the mouse brain. We sequenced 96 retrogradely labeled cells and 295 transsynaptically labeled cells using single-cell RNA-seq, and 4130 retrogradely labeled cells and 2914 transsynaptically labeled cells in situ. We found that the transcriptomic identities of rabies virus-infected cells can be robustly identified using both single-cell RNA-seq and in situ sequencing. By associating gene expression with connectivity inferred from barcode sequencing, we distinguished long-range projecting cortical cell types from multiple cortical areas and identified cell types with converging or diverging synaptic connectivity. Combining in situ sequencing with barcoded rabies virus complements existing sequencing-based neuroanatomical techniques and provides a potential path for mapping synaptic connectivity of neuronal types at scale.

In the brain, messages are relayed from one cell to the next through intricate networks of axons and dendrites that physically interact at junctions known as synapses. Mapping out this synaptic connectivity ­ that is, exactly which neurons are connected via synapses ­ remains a major challenge. Monosynaptic tracing is a powerful approach that allows neuroscientists to explore neural networks by harnessing viruses which spread between neurons via synapses, in particular the rabies virus. This pathogen travels exclusively from 'postsynaptic' to 'presynaptic' neurons ­ from the cell that receives a message at a synapse, back to the one that sends it. A modified variant of the rabies virus can therefore be used to reveal the presynaptic cells connecting to a population of neurons in which it has been originally introduced. However, this method does not allow scientists to identify the exact postsynaptic neuron that each presynaptic cell is connected to. One way to bypass this issue is to combine monosynaptic tracing with RNA barcoding to create distinct versions of the modified rabies virus, which are then introduced into separate populations of neurons. Tracking the spread of each version allows neuroscientists to spot exactly which presynaptic cells signal to each postsynaptic neuron. So far, this approach has been used to examine synaptic connectivity in neurons grown in the laboratory, but it remains difficult to apply it to neurons in the brain. In response, Zhang, Jin et al. aimed to demonstrate how monosynaptic tracing that relies on barcoded rabies viruses could be used to dissect neural networks in the mouse brain. First, they confirmed that it was possible to accurately detect which version of the virus had spread to presynaptic neurons using both in situ and single-cell RNA sequencing. Next, they described how this information could be analysed to build models of potential neural networks, and what type of additional experiments are required for this work. Finally, they used the approach to identify neurons that tend to connect to the same postsynaptic cells and then investigated what these have in common, showing how the technique enables a finer understanding of neural circuits. Overall, the work by Zhang, Jin et al. provides a comprehensive review of the requirements and limitations associated with monosynaptic tracing experiments based on barcoded rabies viruses, as well as how the approach could be optimized in the future. This information will be of broad interest to scientists interested in mapping neural networks in the brain.

Multisensory gamma stimulation promotes glymphatic clearance of amyloid.

The glymphatic movement of fluid through the brain removes metabolic waste1-4. Noninvasive 40 Hz stimulation promotes 40 Hz neural activity in multiple brain regions and attenuates pathology in mouse models of Alzheimer's disease5-8. Here we show that multisensory gamma stimulation promotes the influx of cerebrospinal fluid and the efflux of interstitial fluid in the cortex of the 5XFAD mouse model of Alzheimer's disease. Influx of cerebrospinal fluid was associated with increased aquaporin-4 polarization along astrocytic endfeet and dilated meningeal lymphatic vessels. Inhibiting glymphatic clearance abolished the removal of amyloid by multisensory 40 Hz stimulation. Using chemogenetic manipulation and a genetically encoded sensor for neuropeptide signalling, we found that vasoactive intestinal peptide interneurons facilitate glymphatic clearance by regulating arterial pulsatility. Our findings establish novel mechanisms that recruit the glymphatic system to remove brain amyloid.

Antigen presenting cells in cancer immunity and mediation of immune checkpoint blockade.

Antigen-presenting cells (APCs) are pivotal mediators of immune responses. Their role has increasingly been spotlighted in the realm of cancer immunology, particularly as our understanding of immunotherapy continues to evolve and improve. There is growing evidence that these cells play a non-trivial role in cancer immunity and have roles dependent on surface markers, growth factors, transcription factors, and their surrounding environment. The main dendritic cell (DC) subsets found in cancer are conventional DCs (cDC1 and cDC2), monocyte-derived DCs (moDC), plasmacytoid DCs (pDC), and mature and regulatory DCs (mregDC). The notable subsets of monocytes and macrophages include classical and non-classical monocytes, macrophages, which demonstrate a continuum from a pro-inflammatory (M1) phenotype to an anti-inflammatory (M2) phenotype, and tumor-associated macrophages (TAMs). Despite their classification in the same cell type, each subset may take on an immune-activating or immunosuppressive phenotype, shaped by factors in the tumor microenvironment (TME). In this review, we introduce the role of DCs, monocytes, and macrophages and recent studies investigating them in the cancer immunity context. Additionally, we review how certain characteristics such as abundance, surface markers, and indirect or direct signaling pathways of DCs and macrophages may influence tumor response to immune checkpoint blockade (ICB) therapy. We also highlight existing knowledge gaps regarding the precise contributions of different myeloid cell subsets in influencing the response to ICB therapy. These findings provide a summary of our current understanding of myeloid cells in mediating cancer immunity and ICB and offer insight into alternative or combination therapies that may enhance the success of ICB in cancers.

FinaleMe: Predicting DNA methylation by the fragmentation patterns of plasma cell-free DNA.

Analysis of DNA methylation in cell-free DNA (cfDNA) reveals clinically relevant biomarkers but requires specialized protocols and sufficient input material that limits its applicability. Millions of cfDNA samples have been profiled by genomic sequencing. To maximize the gene regulation information from the existing dataset, we developed FinaleMe, a non-homogeneous Hidden Markov Model (HMM), to predict DNA methylation of cfDNA and, therefore, tissues-of-origin directly from plasma whole-genome sequencing (WGS). We validated the performance with 80 pairs of deep and shallow-coverage WGS and whole-genome bisulfite sequencing (WGBS) data.

Leveraging single-cell ATAC-seq and RNA-seq to identify disease-critical fetal and adult brain cell types.

Prioritizing disease-critical cell types by integrating genome-wide association studies (GWAS) with functional data is a fundamental goal. Single-cell chromatin accessibility (scATAC-seq) and gene expression (scRNA-seq) have characterized cell types at high resolution, and studies integrating GWAS with scRNA-seq have shown promise, but studies integrating GWAS with scATAC-seq have been limited. Here, we identify disease-critical fetal and adult brain cell types by integrating GWAS summary statistics from 28 brain-related diseases/traits (average N = 298 K) with 3.2 million scATAC-seq and scRNA-seq profiles from 83 cell types. We identified disease-critical fetal (respectively adult) brain cell types for 22 (respectively 23) of 28 traits using scATAC-seq, and for 8 (respectively 17) of 28 traits using scRNA-seq. Significant scATAC-seq enrichments included fetal photoreceptor cells for major depressive disorder, fetal ganglion cells for BMI, fetal astrocytes for ADHD, and adult VGLUT2 excitatory neurons for schizophrenia. Our findings improve our understanding of brain-related diseases/traits and inform future analyses.

Integrating single-cell RNA-seq datasets with substantial batch effects.

Integration of single-cell RNA-sequencing (scRNA-seq) datasets has become a standard part of the analysis, with conditional variational autoencoders (cVAE) being among the most popular approaches. Increasingly, researchers are asking to map cells across challenging cases such as cross-organs, species, or organoids and primary tissue, as well as different scRNA-seq protocols, including single-cell and single-nuclei. Current computational methods struggle to harmonize datasets with such substantial differences, driven by technical or biological variation. Here, we propose to address these challenges for the popular cVAE-based approaches by introducing and comparing a series of regularization constraints. The two commonly used strategies for increasing batch correction in cVAEs, that is Kullback-Leibler divergence (KL) regularization strength tuning and adversarial learning, suffer from substantial loss of biological information. Therefore, we adapt, implement, and assess alternative regularization strategies for cVAEs and investigate how they improve batch effect removal or better preserve biological variation, enabling us to propose an optimal cVAE-based integration strategy for complex systems. We show that using a VampPrior instead of the commonly used Gaussian prior not only improves the preservation of biological variation but also unexpectedly batch correction. Moreover, we show that our implementation of cycle-consistency loss leads to significantly better biological preservation than adversarial learning implemented in the previously proposed GLUE model. Additionally, we do not recommend relying only on the KL regularization strength tuning for increasing batch correction, as it removes both biological and batch information without discriminating between the two. Based on our findings, we propose a new model that combines VampPrior and cycle-consistency loss. We show that using it for datasets with substantial batch effects improves downstream interpretation of cell states and biological conditions. To ease the use of the newly proposed model, we make it available in the scvi-tools package as an external model named sysVI. Moreover, in the future, these regularization techniques could be added to other established cVAE-based models to improve the integration of datasets with substantial batch effects.

Dissecting the tumor microenvironment in response to immune checkpoint inhibitors via single-cell and spatial transcriptomics.

Cancer is a disease that undergoes selective pressure to evolve during its progression, becoming increasingly heterogeneous. Tumoral heterogeneity can dictate therapeutic response. Transcriptomics can be used to uncover complexities in cancer and reveal phenotypic heterogeneity that affects disease response. This is especially pertinent in the immune microenvironment, which contains diverse populations of immune cells, and whose dynamic properties influence disease response. The recent development of immunotherapies has revolutionized cancer therapy, with response rates of up to 50% within certain cancers. However, despite advances in immune checkpoint blockade specifically, there remains a significant population of non-responders to these treatments. Transcriptomics can be used to profile immune and other cell populations following immune-checkpoint inhibitor (ICI) treatment, generate predictive biomarkers of resistance or response, assess immune effector function, and identify potential immune checkpoints. Single-cell RNA sequencing has offered insight into mRNA expression within the complex and heterogeneous tumor microenvironment at single-cell resolution. Spatial transcriptomics has enabled measurement of mRNA expression while adding locational context. Here, we review single-cell sequencing and spatial transcriptomic research investigating ICI response within a variety of cancer microenvironments.

Structural variants involved in high-altitude adaptation detected using single-molecule long-read sequencing.

Structural variants (SVs), accounting for a larger fraction of the genome than SNPs/InDels, are an important pool of genetic variation, enabling environmental adaptations. Here, we perform long-read sequencing data of 320 Tibetan and Han samples and show that SVs are highly involved in high-altitude adaptation. We expand the landscape of global SVs, apply robust models of selection and population differentiation combining SVs, SNPs and InDels, and use epigenomic analyses to predict enhancers, target genes and biological functions. We reveal diverse Tibetan-specific SVs affecting the regulatory circuitry of biological functions, including the hypoxia response, energy metabolism and pulmonary function. We find a Tibetan-specific deletion disrupts a super-enhancer and downregulates EPAS1 using enhancer reporter, cellular knock-out and DNA pull-down assays. Our study expands the global SV landscape, reveals the role of gene-regulatory circuitry rewiring in human adaptation, and illustrates the diverse functional roles of SVs in human biology.

An encyclopedia of enhancer-gene regulatory interactions in the human genome.

Identifying transcriptional enhancers and their target genes is essential for understanding gene regulation and the impact of human genetic variation on disease1-6. Here we create and evaluate a resource of >13 million enhancer-gene regulatory interactions across 352 cell types and tissues, by integrating predictive models, measurements of chromatin state and 3D contacts, and largescale genetic perturbations generated by the ENCODE Consortium7. We first create a systematic benchmarking pipeline to compare predictive models, assembling a dataset of 10,411 elementgene pairs measured in CRISPR perturbation experiments, >30,000 fine-mapped eQTLs, and 569 fine-mapped GWAS variants linked to a likely causal gene. Using this framework, we develop a new predictive model, ENCODE-rE2G, that achieves state-of-the-art performance across multiple prediction tasks, demonstrating a strategy involving iterative perturbations and supervised machine learning to build increasingly accurate predictive models of enhancer regulation. Using the ENCODE-rE2G model, we build an encyclopedia of enhancer-gene regulatory interactions in the human genome, which reveals global properties of enhancer networks, identifies differences in the functions of genes that have more or less complex regulatory landscapes, and improves analyses to link noncoding variants to target genes and cell types for common, complex diseases. By interpreting the model, we find evidence that, beyond enhancer activity and 3D enhancer-promoter contacts, additional features guide enhancerpromoter communication including promoter class and enhancer-enhancer synergy. Altogether, these genome-wide maps of enhancer-gene regulatory interactions, benchmarking software, predictive models, and insights about enhancer function provide a valuable resource for future studies of gene regulation and human genetics.

Power of inclusion: Enhancing polygenic prediction with admixed individuals.

Admixed individuals offer unique opportunities for addressing limited transferability in polygenic scores (PGSs), given the substantial trans-ancestry genetic correlation in many complex traits. However, they are rarely considered in PGS training, given the challenges in representing ancestry-matched linkage-disequilibrium reference panels for admixed individuals. Here we present inclusive PGS (iPGS), which captures ancestry-shared genetic effects by finding the exact solution for penalized regression on individual-level data and is thus naturally applicable to admixed individuals. We validate our approach in a simulation study across 33 configurations with varying heritability, polygenicity, and ancestry composition in the training set. When iPGS is applied to n = 237,055 ancestry-diverse individuals in the UK Biobank, it shows the greatest improvements in Africans by 48.9% on average across 60 quantitative traits and up to 50-fold improvements for some traits (neutrophil count, R2 = 0.058) over the baseline model trained on the same number of European individuals. When we allowed iPGS to use n = 284,661 individuals, we observed an average improvement of 60.8% for African, 11.6% for South Asian, 7.3% for non-British White, 4.8% for White British, and 17.8% for the other individuals. We further developed iPGS+refit to jointly model the ancestry-shared and -dependent genetic effects when heterogeneous genetic associations were present. For neutrophil count, for example, iPGS+refit showed the highest predictive performance in the African group (R2 = 0.115), which exceeds the best predictive performance for the White British group (R2 = 0.090 in the iPGS model), even though only 1.49% of individuals used in the iPGS training are of African ancestry. Our results indicate the power of including diverse individuals for developing more equitable PGS models.

Multitissue H3K27ac profiling of GTEx samples links epigenomic variation to disease.

Genetic variants associated with complex traits are primarily noncoding, and their effects on gene-regulatory activity remain largely uncharacterized. To address this, we profile epigenomic variation of histone mark H3K27ac across 387 brain, heart, muscle and lung samples from Genotype-Tissue Expression (GTEx). We annotate 282 k active regulatory elements (AREs) with tissue-specific activity patterns. We identify 2,436 sex-biased AREs and 5,397 genetically influenced AREs associated with 130 k genetic variants (haQTLs) across tissues. We integrate genetic and epigenomic variation to provide mechanistic insights for disease-associated loci from 55 genome-wide association studies (GWAS), by revealing candidate tissues of action, driver SNPs and impacted AREs. Lastly, we build ARE-gene linking scores based on genetics (gLink scores) and demonstrate their unique ability to prioritize SNP-ARE-gene circuits. Overall, our epigenomic datasets, computational integration and mechanistic predictions provide valuable resources and important insights for understanding the molecular basis of human diseases/traits such as schizophrenia.